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Abstract 

We construct a model of spontaneous CP violation in Eg supersymmetric grand 
unified theory. In the model, we employ an SU(2)f flavor symmetry and an anoma- 
lous U(1)a symmetry. The SU(2)f flavor symmetry is introduced to provide the 
origin of hierarchical structures of Yukawa coupling and to ensure the universality 
of sfermion soft masses. The anomalous U(1)a symmetry is introduced to realize 
^ ■ the doublet-triplet mass splitting, to provide the origin of hierarchical structures of 

Yukawa couplings, and to solve the \x problem. In the model, CP is spontaneously 
broken by the SU(2)p breaking in order to provide a Kobayashi-Maskawa phase and 
to evade the supersymmetric CP problem. However, a naive construction of the 
model generally leads to unwanted outcome, Arg[/i&*] = 0(1), when CP violating 
effects in the flavor sector are taken into account. We cure this difficulty by impos- 
ing a discrete symmetry and find that this prescription can play additional roles. It 
ensures that realistic up-quark mass and Cabibbo angle are simultaneously realized 
without cancellation between 0(1) coefficients. Also, severe constraints from the 
chromo electric dipole moment of the quark can be satisfied without destabilizing 
the weak scale. The discrete symmetry reduces the number of free parameters, 
but the model is capable of reproducing quark and lepton mass spectra, mixing 
angles, and a Jarlskog invariant. We obtain characteristic predictions V u b ~ 0(A 4 ) 
(A = 0.22) and iV^&Yfil = \Y C \ at the grand unified theory scale. 
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1 Introduction 

Low-energy supersymmetry (SUSY) is one of the leading candidates for physics beyond 
the standard model (SM) [I]. In particular, the minimal supersymmetric extensions of 
the SM (MSSM) exhibit remarkable coincidence of three SM gauge coupling constants 
near 10 16 GeV, and the supersymmetric grand unified theory (SUSY GUT) appears to 
be plausible model behind the MSSM. In these models, however, if generic soft SUSY 
breaking terms are introduced, these terms can induce flavor changing neutral current 
(FCNC) processes and CP violating observables like electric dipole moments (EDMs) 
that are too large to be consistent with the present experiments [2]. These issues are 
called the SUSY flavor problem and the SUSY CP problem, respectively. In this paper, 
we address the later problem in the context of SUSY GUT. We choose E$ as a unification 
group and introduce an SU(2) flavor symmetry. Based on this framework, we construct 
a model of spontaneous CP violation (SCPV) in the flavor sector to provide the origin of 
the Kobayashi-Maskawa (KM) phase and to evade the SUSY CP problem. 

Among several possible grand unification groups [3J SI E] , the E 6 group with SU(2) F 
flavor symmetry exhibits a numbers of attractive features [3 El [7]. 

1. Three fundamental representations of E 6 , 27!, 27 2 , and 273, contain all of the SM 
fermions and the right-handed neutrinos. 

2. Various realistic hierarchical structures of masses and mixings of quarks and lep- 
tons can be naturally derived from the one basic hierarchy of E$ invariant Yukawa 
couplings. 

3. This basic hierarchy results from the SU(2)p flavor symmetry and its breaking ef- 
fects0 

4. The SUSY flavor problem can be avoided by the SU(2)p flavor symmetry without 
destabilizing the weak scale@ To be more concrete, at the leading order, the model 
has the following sfermion soft mass matrices at the GUT scale, as a consequence 
of SU(2)p symmetry [8]. 

\ / ml \ 

ml o , ml = ml (1.1) 
m 2 3 J \ ml ) 

(Here, m 10 and mg are the mass matrices of the 10 and 5 fields of SU(5), where 
the SM matter fields are embedded in these multiplets in a standard way [3].) In 
this paper, we call this form of sfermion soft mass matrix "modified universality" [8j 
EJ [101 [HI [12]. In this set up, if the gluino mass (M 3 ), the up- type Higgs mass 
{ m Hu)i and the supersymmetric Higgs mass (//) as well as m 3 are the weak scale, 

* x An anomalous U(1)a symmetry which we discuss in Sec. 3 also plays a role to generate hierarchical 
structures of Yukawa couplings. 

* 2 Here we assume that the D-term contributions are small. A possible way to justify this assumption is 
by adopting non-Abelian discrete symmetry instead of the gauged symmetry which we use in this paper. 
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then the weak scale can be stabilized^ Also, even if there are some deviations 
from the "modified sfermion universality" in flavor space, the FCNC constraints 
can be satisfied by raising m , and this does not destabilize the weak scale as long 
as m 3 is maintained in the weak scale0 Note that this prescription is more natural 
compared to the case for the deviation from the "universality" that is assumed in 
minimal supergravity where is set equal to m®. 

Now, let us move on to the SUSY CP problem. In order to avoid the SUSY CP 
problem, it is desirable that the dimensionful soft SUSY breaking parameters and the 
/i parameter are all real. On the other hand, the Yukawa couplings should be complex 
parameters so as to provide a KM phase [16]. Is it possible to realize these observations 
naturally? To answer this question, in this paper we examine the idea of spontaneous CP 
violation (SCPV) PH EE] within the context of E 6 xSU(2) F SUSY GUT. First of all, by 
imposing CP symmetry, all the coupling constants of the original Lagrangian, including 
soft SUSY breaking terms, can be taken to be real. Then we make a model in which CP is 
spontaneously broken by complex vacuum expectation values (VEV) of flavon fields; i.e., 
the SU(2)p flavor symmetry breaking is also responsible for the CP violation. Here, once 
renormalizability is assumed, it seems difficult to generate complex Yukawa couplings. 
Therefore we allow higher dimensional operators so that complex dimensionless (Yukawa) 
couplings can be induced via nonrenormalizable interaction terms. Since CP is broken 
in the SU(2)f sector, we can expect that the SU(2)p singlet SUSY breaking parameters 
can remain real after the SCPV, contrary to the Yukawa couplings. Also, even if flavor 
nonsinglet SUSY breaking parameters acquire CP phases, their effects on the CP violating 
observables can be suppressed in the current framework by taking m much larger than 
the weak scaleQ Therefore, SCPV provides a way to solve the SUSY CP problem. 

In the previous papers of Ref. [19], an anomalous U(1)a gauge symmetry [201 [2TJ was 
introduced to solve the doublet-triplet splitting (DTS) problem and to provide the origin 
of the hierarchal structures of Yukawa couplings. The U(1)a symmetry was also applied 
to solve the \x problem [22] • However, as we will see in the following section, if we naively 
apply the U(1)a symmetry to the model of SCPV, it generally leads to an unwanted out- 
come, Arg[fib*} = 0(1), and gives rise to the SUSY CP problem. We cure this difficulty 
by imposing discrete symmetry. Interestingly, this discrete symmetry can play additional 
roles in the model. It ensures that realistic up-quark mass and Cabibbo angle are si- 
multaneously realized without cancellation between 0(1) coefficients. Also, the severe 
constraints from the chromo-EDMs (CEDM) can be satisfied without destabilizing the 
weak scale. As a result of the discrete symmetry, the number of free parameters of the 
model is reduced. However, the model is capable of reproducing quark and lepton mass 

* 3 In order to satisfy the lightest CP-even Higgs boson mass bound set by LEPI, one may cm- 
ploy maximal mixing of the stop sector |13j . Another interesting option is light Higgs scenario or 
inverted hierarchy where one assumes that the heaviest CP-even Higgs boson behaves as the SM Higgs 
boson [14] . 

* 4 However, the muon g — 2 anomaly is not explained if tuq >lTeV, since SUSY contributions decouple 

5 If 77i3 is placed in the weak scale in view of the naturalness, there are nondecoupling SUSY contri- 
butions to the CP violating observables. These effects can be suppressed if CP is not effectively broken 
in the up-type (s)quark sector |llj . We will comment this point in Sees. 4 and 5. 
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spectra, mixing angles, and a Jarlskog invariant. We also obtain characteristic predictions 
V ub ~ C(A 4 ) (A = 0.22) and \Y b V cb \ = \Y C \ at the GUT scale. 

This paper is organized as follows. In the next section, we give a brief review of the 
model of E 6 SUSY GUT and SU(2) F flavor symmetry. In Sec. 3, an anomalous U(1)a 
gauge symmetry and its specific vacua are summarized. Based on these two sections, we 
show a simple model of the SCPV that is caused by the VEVs of flavon fields, in Sec. 4. In 
the same section, we also discuss CP violation effects on /x- and 6-terms and indicate that 
they lead to Arg [/icV] = 0(1). Then, we introduce discrete symmetry in order to solve 
this difficulty In Sec. 5, we summarize a model of SCPV in E 6 SUSY GUT with SU(2) F , 
U(1)a, and discrete symmetries. Then, we examine the consequent Yukawa couplings 
and derive some predictions of the model. The last section is devoted to a summary and 
discussion. 



2 E 6 xSU(2) F SUSY GUT 

In this section we briefly summarize the model of E 6 xSU(2) F SUSY GUT discussed in 
Refs. [7J IB]- We mainly focus on the Yukawa sector and sfermion soft masses of the model. 

2.1 E 6 GUT 

In this subsection, we discuss the Yukawa sector of E 6 SUSY GUT. Here the twisting 
mechanism among SU(5) 5 fields plays important roles [SI El IE]- 

For the E 6 group, 27 is the fundamental representation, and in terms of E 6 DSO(IO) xU(l)y 
(and [SO(10)dSU(5)xU(1) v ]) it is decomposed afl 

27 = 16i[10_i + 5 3 + 1_ B ] + 10_ 2 [5 2 + 5'_ 2 ] + U[1' ] . (2.1) 

As shown in (12.11) . each 27 incorporates two 5's(and l's) of SU(5). This nature allows 
a model to simply produce different and realistic Yukawa structures of quark and lepton 
from a single hierarchical structure of an E 6 invariant Yukawa coupling [SI E] • 
Let us introduce the following superpotential: 

W Ee> D y;fvl/f vl/J 7 // 27 + Y?y 2r y 27 C 27 (2.2) 

and assume that the original Yukawa hierarchies are 

A 6 A 5 A 3 



Yg - Y% ~ [ A 5 A 4 A 2 I . (2.3) 
A 3 A 2 1 



Here \I/ 2T (i — 1 — 3) are matter fields, and H 27 and C 27 are Higgs fields that break E 6 
into SO(10) and SO(10) into SU(5), respectivelyH We assume that the MSSM Hide's 



* 6 Here acutes are used to distinguish different 5(l)s. Each subscript represents corresponding U(l) 
charge. 

* 7 In the following discussion, we also introduce antifundamental fields, H 27 and C 27 , to maintain cor- 
responding D-flatness conditions. Also, the adjoint field A 78 is introduced to realize the DTS, employing 
the Dimopoulos-Wilczek (DW) mechanism [53]. A complete list of field contents is given in SccfS] 
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doublets are included in H and C . Hereafter, we parametrize the several hierarchical 
structures of couplings and VEVs of various fields in terms of the Cabibbo angle A. (We 
set A = 0.22.) In Eq. fl2.3p . we only parametrize the hierarchical structures, and the 0(1) 
coefficients are suppressed for the moment. 

In Eq. fl2.ip . once H 27 and C 27 acquire VEVs in the components of l^l'o) and 16i(l_ 5 ), 
respectively, a superheavy mass matrix of rank 3 among 5j, 5[, and 5j of ty 27 is induced, 
through the Yukawa coupling of Eq. (I2.2[) . Therefore the 3 degrees of freedom among 
5' and 5; decouple at the GUT scale. In consequence, the up-type quark Yukawa cou- 
pling (Yfj) remains in the original form of Eq. fj2.3j) but the down-type quark (Yd) and 
charged lepton Yukawa couplings (Ye) differ from it. Importantly, the three massless 
modes among 5^ and 5$ mainly originate from the first two generations of \I/ 27 because the 
third generation 5' 3 and 83 from ^I/ 27 have large Yukawa couplings, i.e., large GUT scale 
masses. This feature is also important for making sfermion universality for the 5 sector, 
as discussed in the following subsection. As an example, when (C)/{H) ~ A ' 5 , we end 
up with the following milder hierarchies for Ye and Ye [7]: 



Note that the hierarchies of Eqs.f l2.3p and f !2.4p are adequate to reproduce several mass 
spectra and flavor mixing angles of quarks and leptonsO However, an exception is the 
up-quark Yukawa coupling Y u . The hierarchy of Eq. ( I2.3[) leads to Y u ~ A 6 but the realistic 
value is Y u ~ A 8 Q We will get back to this issue in Sec 14.31 

2.2 SU(2) flavor symmetry and modified universality 

In this subsection, we briefly discuss the effect of SU(2)p flavor symmetry on the sfermion 
soft masses and see the emergence of the modified universal form of Eq. fll.il) [8]. 

First of all, we assume that the original hierarchy of Eq. (12 .3p (partly) originates from 
the flavor symmetry breaking effects. In order to naturally incorporate the 0(1) top 
quark Yukawa coupling, here, we employ SU(2)p flavor symmetry and treat the first two 
generations of matter fields as the doublet ty 27 whereas the third generation field ^> 2 and 
all the Higgs fields are treated as singlets under SU(2)p. We also introduce flavon fields 
F a and F a that are singlets under Eq and a doublet and an anti-doublet under SU(2)p, 
respectively. Then, as shown later, F a and F a acquire VEVs and break the the SU(2)p. 
These effects generate hierarchical structures of the Yukawa coupling of Eq. (12.31) . 

Suppose that the soft SUSY breaking terms are mediated to the visible sector above 
the scale where Eq and SU(2)p symmetries are respected, such as in gravity mediation. 
Then, ty 27 acquires a soft mass (mo) that is different from the mass of ^> 27 (777,3), m 
general, in an E 6 xSU(2)f symmetric way. As a result, this guarantees whole sfermion 
soft mass degeneracy at the GUT scale, except for 10 3 , as shown in E q.f|l.lD . because 
10i, 10 2 , and all the MSSM 5 fields originate from an identical field frg 7 P°T 

* 8 See Ref.[7] for details. 
* 9 See, e.g., Ref.[lfl[25]. 

* 10 Since possible deviations from Eq. (|2.2[) . which are induced by the SU(2)p and E 6 breaking effects, 




A 6 A 5 ' 5 A 5 
A 5 A 4 - 5 A 4 
A 3 A 2 ' 5 A 2 




(2.4) 
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3 Anomalous U(l) gauge symmetry 



In the previous papers of Ref . [T9] . an anomalous U(1)a gauge symmetry [20], [2T] was 
introduced to solve the doublet-triplet mass splitting and to provide the origin of the 
hierarchical structures of Yukawa couplings. In this section, we summarize essential points 
and comment on what kinds of fields acquire VEVs, and how the magnitudes of these 
VEVs are determined in the framework of U(1)a- 

We introduce an anomalous U(1)a gauge symmetry whose anomalies are canceled by 
the Green-Schwarz mechanism [26J. The theory possesses the Fayet-Iliopoulos term (£ 2 ), 
and we assume its magnitude as £ = AA. (Here A is the cutoff scale of the theory and 
we set A = 1.) Let us denote the symmetries of the theory, except for U(1)a, as G a . 
Then it is shown in Ref. [19] that a system consisting of all the G a and U(1)a invariant 
termJ*"1 has a supersymmetric vacuum where all the fields that are negatively charged 
under U(1)a get VEVs in the following way! 12 



Here Zj denotes the G a singlet field, but the argument can be extended to the case where 
Zi are composite operators that are made by G a nonsinglet fields. For example, in case the 
where Z~ = XX (here X denotes the G a nonsinglet field and X is the antirepresentation 
of X), (XX) = leads to (X) = (X) = A"^^ 2 , once the D-flatness condition of 

G a is taken into account. 

A U(1)a symmetry and its specific SUSY vacuum are applied in several aspects of 
phenomenological model building. For example, an appropriate U(1)a charge assignment 
for the Higgs sector can ensure the DTS via the DW mechanism [191 E3]. Here note that 
since a positively charged field does not acquire a VEV in this U(1)a vacuum, only the 
terms that are linear in a positively charged field are relevant to the determination of the 
symmetry breaking structure of the model. 

It is also important to mention that a term whose total U(1)a charges are negative 
does not appear at the U(1)a breaking vacuum. The reason is that this type of term 
should originally be accompanied by at least one positively charged field but its VEV 
is always vanishing according to Eq. fl3.ip (SUSY- zero mechanism) [271 US]- Importantly, 
this property can be applied to solve the \x problem [22]. To be more precise, when the 
U(1)a charges of the Higgs fields are set to be negative, the supersymmetric /i-term is 
forbidden by the SUSY-zero mechanism. However, this SUSY-zero mechanism will be 
broken by the amount of the SUSY breaking scale, once soft SUSY breaking effects are 
taken into account in the total Lagrangian. So the appropriate scale of the //-term can 
be induced. Employing these properties of U(1)a, in the following sections we construct 



are argued in Ref. [8], we do not repeat the detailed argument here but the FCNC constraints can be 
satisfied in the current model by employing the argument given in the Introduction. 

* n Here we include all the nonrenormalizable operators. Also, each term is assumed to have an 0(1) 
coefficient. 

* 12 From now on, each superfield is denoted by an uppercase letter, whereas the corresponding lower 
case letter indicates an associated U(1)a charge. The consistency of Eq. (|3.f[) requires the number of 
positively charged fields to be larger than that of negatively charged fields by one. 




and examine a model of SCPV in the Er SUSY GUT. 
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4 Spontaneous CP violation 



In this section, we first discuss the mechanism of SCPV that is caused by the SU(2)p 
breaking sector. After specifying a source for the CP violating phase, we examine how 
this phase will appear in the parameters of the MSSM. Particularly, we refer to a solution 
for the \x problem in the context of U(1)a and indicate that it leads to an unwanted 
outcome, Arg[/i&*] = 0(1), once CP violating effects are taken into account. Next, we 
discuss discrete symmetry in order to solve this difficulty. 



4.1 Spontaneous CP violation 

In this subsection, adopting a U(1)a symmetry and its characteristic vacuum, we discuss 
the mechanism of SCPV that is caused by the SU(2)p breaking sector. 

Let us introduce the EexSU(2)p singlet field S(s > 0). Then, in general, we obtain 
the following superpotential made of F a (f < —1), F a (f < —1), and S: 

W s = X S S[J2 c n X if+f)n {FF) n ] (4.1) 

n=0 

We assume that 0(1) coefficients c n are all real and that the theory is originally CP 



invariantj* 13 l n are natural numbers which are truncated at rif, where s + (/ + f){ n f + 
1) becomes negative (SUSY-zero mechanism). The appearance of A in Eq. fl4.ip can be 
understood as originating from the VEV of an operator whose U(1)a charge is —1. From 
Eq.f l4.ip . the F- flatness condition with respect to S leads to 

A s [c + Cl A (/+/) (FF) + • • ■ + c n/ A n ' (/+/) (FF) n /] = . (4.3) 

Thus, when nf > 2, it generally forces (FF) to acquire an imaginary phase. Therefore, 
at this stage the CP invariance is spontaneously broken. Using SU(2) F gauge symmetry 
and taking its D-flatness condition into account, hereafter we adopt the following basis 
where only (F2) acquires an imaginary phase: 

(Fa) ~ I eipX -if+H ) > ~ ( A _2+I> I (4-4) 



4.2 Spontaneous CP violation and /i-term generation 

In this subsection, we discuss the effect of Eq fl4.4l) on the generation of the /1- and 6-terms. 
Here an appropriate //-term is induced by the SUSY breaking effect on the SUSY-zero 
mechanism] 14 

* 13 In general, following terms can also appear in the bracket of Eq. (|4.1|) . 

However, these terms are not important for the determination of (FF). From Eq. (|3.1j) . it is understood 
that their signatures are summarized in cq of Eq. (|4.ip . 

14 Note that, within the context of the U(1)a summarized in the previous section, the Giudice-Masiero 
mechanism [25] does not work well since it requires \2h\ < 1 [29] , 
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Let us assume that the MSSM Higgs doublets are contained in 5 and 5' of H 27 . 
Then if the U(1)a charge of H 27 is negative (h<0), there is no supersymmetric mass 
term for the MSSM Higgs doublets. This is because the positively charged field S does 
not acquire VEVs in the SUSY limit (SUSY-zero mechanism). Here, a possible term is 
X s+3h SH 27 H 27 H 27 -»■ A s+3 ' i (5)(i/ 1 )i/ 10 i/ 10 , but (S) is vanishing in the SUSY limit. 

However, the situation changes once we take SUSY breaking effects into account in 
the total Lagrangian [30j 122] . For a schematic explanation, we first discuss a simplified 
version of the system which consists of EqxSU(2)f singlet fields S(s > 0) and Z(z < — 1). 
The superpotential is given as 

W D X s S + \ S+Z SZ + . . . (4.5) 



and we also introduce the corresponding general soft SUSY breaking terma* 15 



V mft = E rh 2 x \X\ 2 + A s X s S + A z X s+z SZ+.... (4.6) 

x=s,z 

Here we assume that all the soft parameters rhs, friz, As and Az are the weak scale. As 

usual, the whole potential is composed of V = Yl Wx\ 2 + V so f t JJj where 

x=s,z 

\F S \ 2 = & 2 = \X S + X S+Z Z\ 2 , \F Z \ 2 = & 2 = \X S+Z S\ 2 . (4.7) 
do oZ 

Here note that, in the SUSY limit, (Z) is determined so that \F$\ of Eq.f l4.7p vanishes. 
Also, (S) = satisfies the vanishing \F Z \. However, once SUSY breaking effects are taken 
into account, there is no need for F-flatness conditions to be satisfied. Actually, the total 
potential V which includes SUSY breaking effects makes (S) and (F$) nonvanishing. First 
of all, the extremum condition of V with respect to S determines (S) as 

(S)~X^ s+2z \A s + A z ). (4.8) 

Here the shift of (S) from zero can be understood as the balance between the SUSY mass 



term and the tadpole term that is induced by the SUSY breaking effect J 17 l Also, F$ is 
fixed by the extremum condition of V with respect to AZ, where AZ is defined as the 
deviation of (Z) from its SUSY VEV as (Z) = X~ z + AZ. 



5V 
SAZ 



F*X S+Z + m\X~ z + A Z X S+Z S = 

(F*) = A" (s+22) m| + A Z (S) (4.9) 



* 15 In Eq.gSJ), abbreviation denotes X 2s+z S 2 Z, X s+2z SZ 2 , etc. Since SUSY breaking effects does not 
largely change Eq(3.1), these terms are not relevant and omitted in the following arguments. This is the 



same for Eq. (|4.6|) . 

* 16 Note that the U(1)a D-flatness is satisfied by the VEV of an operator whose U(1)a charge is — 1. 

* 17 We plot the potential respect to S in FigiTJ Here the dotted parabola shows the SUSY mass term 
(|A S+Z 5| 2 ) and the dashed-dotted line indicates the SUSY breaking tadpole contribution (AsX s S + 
AzX s+z S(Z)). The solid line is the sum of these contributions. 
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Figure 1: SUSY breaking effect on the potential of S 



Combining Eqs.( l4.8j ) and (14. 9p . the /i- and b-terms are induced as follows: 

£ D X s+3h SHH(H)\ e 2 + \ s+3h A SH SHH(H) 

\ s + 2h ^((S)HH\ e2 + ((F s ) + (S)A SH )HH) (4.10) 

Therefore, when z is appropriately chosen, the /i and b parameters appear as the weak 
scale. 

Now, let us examine the effects of the CP violation of Eq. (l4.4p on the generation of 
the n and b parameters. So far we have assumed that all the soft SUSY breaking terms 
and 0(1) coefficients were real. Therefore, the /i and b parameters are real. However, this 
result can be changed, in general, once flavon fields acquire complex VEVs. In order to 
see this, we extend Eq. (l4.5p to include Eq. fl4.ip . Also, we take the corresponding SUSY 
breaking terms 



A\ s S[J2c' n X (f+I)n (FF) n } (4.11) 



n=0 



into account. Here note that each 0(1) coefficient c' n appearing in Eq. (14.111) is different 



from that of Eq. fl4.ll) . in general] 18 l Therefore 



AX s S{c' + ^X f+f {FF) + ■■■ + Z(/ nf \V+V n f{FF) n f} (4.12) 

is generically nonvanishing, even though the left-hand side of (14. 3p vanishes. Note that, 
from Eq. fl3.ip . each higher dimensional operator in Eq. fl4.12l) gives a contribution that 
is of the same order of the original tadpole term, c' Q AX s S, but the effective coefficient is 
complex. Therefore, Eq.( 14.1l|) turns into A$X S S of Eq.( 14.6p . where As is complex and 
Arg[A 5 ] ~ 0(1). Then, we conclude from Eqs.(@~8]), and fl4~T0D that the fi and b 

parameters turn out to be complex parameters and Arg[/i6*] = 0(1), in general. 



4.3 SUSY CP problem and discrete symmetry 

As we saw in the previous subsection, CP violating VEVs of flavon fields generically lead 
to complex /i and b parameters and Arg[/i6*] = 0(1). This rephasing invariant complex 

* 18 When each 0(1) coefficient appearing in a superpotential term equals to the corresponding coefficient 
of a soft SUSY breaking term, the above mechanism for the generation of ^-term does not work. 
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phase induces (C)EDMs of quarks and leptons, and gives rise to the SUSY CP problem. 
Though the decoupling procedure can reduce the SUSY contribution to the (C)EDMs, 
we discuss another way to avoid this issue in order not to destabilize the weak scale. 

As we discussed below Eq. fl4.lip . the couplings among S and flavon fields led to 
complex fi and b parameters and Arg[yu6*] = 0(1). Therefore, we introduce a discrete 



symmetry so that coupling among S and FF is forbidden]* 19 ! Specifically, we exploit 
the discrete symmetry in order to discriminate F from the others assigning a nontrivial 
discrete charge for F, since only F acquires a complex VEV in the basis defined in Eq. (I4.4p . 
One may think that the introduction of an additional symmetry seems ad hoc, but we can 
find the utility of this discrete symmetry from several different point of view, as follows. 

As we mentioned at the end of Sec.2.1, a naive assumption for the Yukawa hierarchy 
of Eq. fl2.3p leads to a relatively large up-quark Yukawa coupling. Therefore, there should 
be cancellation between 0(1) coefficients of Eq. fl2.3p . in order to reproduce realistic up- 
quark Yukawa coupling (Y u ~ A 8 ) [211 [25]. However, the situation can be improved if the 
following types of interactions are only responsible for the generation of Y# : 

Q y2ip a +a^,a27 ^78^27 q \ 

^2^+^27^78^27 ^ a +f) pa^,27 pb^27 yPa+^i+J pa^, ^ ^ | yh jj27 

A^+^ 3+ % 3 F a ^ A 2 ^3^ 3 / 

/ Q B _ L A 5 *i* 2 \ 

-> Q B _ L A 5 ^i A 4 * 2 * 2 A 2 * 2 *3 H (4.13) 

\ A 2 $ 3 ^2 *3*3 / 

Here A 78 is assumed to acquire B — L conserving VEVs [19] and Qb~l represents the B — L 
charge of the corresponding component of \1/ 27 . Then, Eq. fl4.13p leads to Y u ~ Qb-l^ 
which can be an order of magnitude smaller than the original expectation, Y u ~ A 



20 



Simply, Eq. fl4.13p is ensured if there is a discrete symmetry which discriminates F from 



■21 



the others and forbids F a ^ 2 J F b ^ 27 H 27 and F a mf F b ^ 27 H 27 

Also, another utility of discrete symmetry arise from the examination of the SUSY CP 
problem that is caused even if the \x- and tVterms are real. In the case where the modified 
sfermion universality of Eq.f ll.ip is adopted, the nonuniversality of the up-type squark 
sector (m^o) induces quark a CEDM [2U[TT]. This CEDM exceeds experimental bound if 
m 3 is placed in the weak scale. Here, the CP violating phase is provided by the up-type 
quark Yukawa coupling, even if /x, b, M3, and A parameters are real. Therefore, when we 
respect the weak scale stability, the up-type quark Yukawa coupling should also be real 



* 19 In actual model building we should introduce an additional positive singlet field S" which plays a 
role of S appeared in Eq. (|4.1[) . Then this S' may acquire complex VEV in its scalar and F components. 
However, the coupling between S' and H 3 can be forbidden by an appropriate symmetry. Or, since the 
magnitude of (S") is dictated by the smallest negative charge of the relevant field, another Z' field can 
make (5") sufficiently small, and ensures y,- and &-terms to be almost real. 

20 In general Yjj and Yjj are the combinations of Yh and Yq. However, we make Higgs sector so that 
the MSSM up-type Higgs doublet Hjj originates solely from H u [H b ] but the down-type Higgs doublet 
Hd originates as Ho = Hd[H 5 ]+A°- 5 i[C 5 ]. (See Appendix \M) In this case, Yjj is combinations of Yh 
and Yc but Yjj is directly given by Yh- 

* 21 Note that, in this case, F a ^ a ^ 3 ->■ A 3 *i* 3 is also forbidden. 
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parameters other than M3, /i, and b \* 22 \ This requirement can also be ensured if there is 
a discrete symmetry which discriminates F from the others and F does not contribute to 
the generation of Yjj, since only F acquires complex VEV, 



*23 



In connection with this type of SUSY CP problem that is caused by the flavor off- 
diagonal complex soft term, the SU(2)p breaking effect can also induce similar flavor 
off-diagonal complex soft term. For example, a higher dimensional term in the Kahler po- 
tential, m 2 (ty 27 y e ab F a ty 27 , induces a complex (m 2 )i3 component of the order of e ip A 3 m 2 , 
and it contributes to the up-quark (C)EDM. However, this effect can also be forbidden if 
there is a discrete symmetry which discriminates F from the others. 

Following these observations, we introduce a discrete symmetry to our model so that 
all the above requirements can be fulfilled. When we adopt a cyclic symmetry, Z3 is the 
smallest group that can allow F a ^ 2 a 7 F b ^ 27 H 27 but forbid F a ^l 7 F b ^> 27 H 27 . Then, e.g., 
F is assigned Z3 = 1 and all the other fields are assigned Z3 = 0, in order to discriminate 
F from the others. Here, an exception is C 27 . As explained in the next section, the model 
can not be realistic without F a ^l 7 F b ^! 27 C 27 and F a ^l 7 ty 27 C 27 , since it supplies the KM 
phase through the down-type quark Yukawa coupling. Therefore, C 27 should be assigned 
Z3 = 2. Once a Higgs field is charged under this new Z3 symmetry, some modifications 
are necessary in the Higgs sector that is originally composed to ensure the DTS [19]. In 
the following discussion, therefore, we extend Z 3 into Z 6 where Z Q incorporates the Z 2 
symmetry that was used to ensure that the adjoint Higgs field acquired B-L conserving 
VEV [321 [TH]. Here, a simple choice of Z 6 charges is Z 6 [F] = 2 and Z 6 [C] = 4. However, 
we set Zq [F] = 1 and Z e [C] = 5 in order to utilize the degrees of freedom of the enlarged 
discrete symmetry. 



5 A model of SCPV in E 6 SUSY GUT 

In this section, we summarize a model of SCPV in the Eq SUSY GUT with SU(2)p, U(1)a, 
and Zq symmetries. After giving brief explanations for the field contents, we focus on the 
effects of the discrete symmetry on the Yukawa couplings of the model. We examine 
consequent Yukawa couplings and derive some predictions from them. 

5.1 Field content 

Here we summarize the field content of the model and its representations under the E 6 , 
SU(2) F , U(1)a and Z e symmetries. 

We introduce the following fields, which are listed in Table [TJ First of all, \I/ 27 is 
a matter field, and we make its first two generations as a doublet (\I/ 27 ) and the third 
generation as a singlet (^f 7 ) of the SU(2)p, respectively. The SU(2)p and CP are simulta- 
neously broken by the VEVs of flavon fields F a and F a as in Eq.f l4.4p . All the other fields 
are singlets under SU(2)p. H 27 is the field whose VEV breaks E 6 into SO(10), and H 27 

* 22 We assume that each Yukawa coupling has a corresponding interaction with the SUSY breaking 
spurion field (X = 6 2 msuSY), like ^i^jH(C), where cf^ C ^ are the real 0(1) coefficient. 

Then if Yf is real (complex), the corresponding Af term is also real (complex). 

* 23 See EH 
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$27 

a 




F a 


pa 


H 27 


^27 


C 27 


£27 


(Ji27 


(JI27 


A 78 


A' 78 




27 


27 


1 


1 


27 


27 


27 


27 


27 


27 


78 


78 


SU{2) F 


2 


1 


2 


2 


1 


1 


1 


1 


1 


1 


1 


1 


U(1)a 


4 


3 
2 


3 
2 


5 
2 


-3 


2 


-4 





7 


9 


-1 


4 


Ze 








1 











5 





3 


3 


3 


3 



^0 ^3 5 
Eg 1111 

St/(2) F 1111 
U(1) A -1 -2 -5 9 
Z 6 3 4 



Table 1: Field contents and charge assignment under E6xSU(2)fxU(1)axZ6 



is introduced to maintain the D-flatness condition. C 27 is the field whose VEV breaks 
SO(10) into SU(5), and C 27 is introduced to maintain the corresponding D-flatness con- 
dition. A 78 and A' 78 are the adjoint fields. Here, the F-flatness conditions with respect 
to A' 78 make A 78 acquire DW-type VEVs for the SO(10) adjoint component to solve the 
DTS problem [23]. C' 27 , C' 27 , and Z are introduced to give masses for the NG modes 
and to maintain DW-type VEVs [32] . 

In Table [U U(1)a charges are assigned so that the DTS and appropriate Yukawa hi- 
erarchies are ensured. Also, Z 6 charges are determined so that the requirement discussed 
in Sec J4.3l can be fulfilled. A more detailed discussion for the symmetry breaking in this 
model can be found in [15] . As a result of the SUSY-zero mechanism and Z 6 discrete 
symmetry, possible forms of Yukawa couplings are characteristically restricted. In the fol- 
lowing subsections, we examine consequent Yukawa couplings and derive some predictions 
of the model. 



5.2 Yukawa couplings 

In this subsection, we examine Yukawa couplings of the model restoring 0(1) coefficients. 

Under the charge assignment of Table [fl the following interactions between matter 
and Higgs fields are allowed. 

dV a (A,Z 3 )V a 

Y H : dty a (A,Z 3 )ty a c\ 2 ^ + ^F a m a F h m b bX^ + ^ + ^F a ^ a ^ 3 ] \ h H (5.1) 
b\^ + ^ + f^ 3 F a ^a aA 2 ^3^ 3 



f\ 2 ^+f+fF a y a F b V b g\^ + ^+fF a ty a ty 3 

Y c : | f\ 2 ^ + f+fF a m a F h m b ] X C C 

g\^+^+f^ 3 F a y a 

(5.2) 
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5' 2 


B's 


5i 


5 2 


5 3 





ad 5 X 5 








/e i5 A 5 ' 5 


ge i5 X 


ad 5 \ 5 


cA 4 


6A 2 


/e i<5 A 5 - 5 











b\ 2 


a 


#e i<5 A 3 - 5 









These interactions are responsible for the generation of Y# and Ye- Here we explicitly 
write in 0(1) coefficients, a, b, c, d, g and /, for later discussion. In Eq.f l5.ip . ^/ a (A, Z 3 )^/ a H 
includes X 2 ^ +a+Z3+h ^ a Z 3 A^ a H , A 2 ^ +a )+ h ^<M 2 # a F, etc. Since A acquires B-L con- 
serving VEVs ((A) ~ Qb-lX), the effects of <i^ a ((A, Z 3 ))^ a H differ for the different 
components of \E^ 7 ', as we reparametrize below. Note that all the 0(1) coefficients are 
assumed to be real numbers because of the original CP symmetry. 

When the Higgs fields and flavon fields acquire VEVs, Eqs. fl5.ll) and (15.21) induce the 
following mass matrix for 5j, 5^ and 5^: 



(H). (5.3) 



Here each power of A is determined by the corresponding U(1)a charge. It is important 
to note that a = 1 for the colored Higgs components (H c , H ) of 5 and 5', and a = for 
the doublet Higgs components (H u ,Hd) of 5 and 5', since the (1,2) and (2,1) elements of 
Eq. (15.31) originate from the B-L conserving VEV of A I* 24 ! Then, we take the following con- 
ventions dV((Z,A))VH D d 5 X 5 e ab H c [^ a \H°[%}l'[H), =^d q X^e ab Q[^ a ]qg[^> b ]H u [H](or 
H d [H\), ^diX^e^alL^^HdiH], where e 12 = -e 21 = 1, e 11 = e 22 = and d 5 ,d q ,d l 
are the 0(1) coefficients that are different from each other in general. In Eq. fl5.3p . since 
a = for H u and Hd components of 5 and 5', the ^[5^] component becomes a purely 
massless mode. In other words, L[5 3 ] does not contain a massless mode whose main mode 
is Hd{5[]. This fact requires the down-type Higgs doublet to be composed of not only 
Hd[H 5 } but also £[C 5 ]; otherwise the determinant of Ye vanishes, as we will see later. 

Now we derive the Yukawa couplings Yu, Ye, and Ye in a semianalytical way, assuming 
that each 0(1) coefficient does not largely alter the hierarchical structures that originate 
in Eqs. fl5.ll) . (15. 2 p and f 15 . 3 j) . First of all, Y v is simply derived from Eq. fl5.ip by extracting 



the corresponding components a 



: 25 







Ui 3 





— -^dqX 





7j dqX^ 


cA 4 


bX 2 





bX 2 


a 



Y v = Q 2 \ H A cA bX ■ ( 5 - 4 ) 



For Yd and Ye, we first derive the relation between the gauge eigen modes and the mass 
eigen modes of 5 from Eq. fl5.3p . Then we replace 5, and 5^ of ty a and \l/ 3 , whic appeared 
in Eqs.f l5.lp and (15. 2p . in terms of massless modes. In the same step, we also replace 
Hd[5' H ] and L[5c] by H D and /3hX°' 5 H d , respectively. Here (3h denotes another 0{\) 
coefficient which enters in this step. When we leave only the leading order contributions, 



* 24 This is true even if higher dimensional term, like Zz>5> a A z ^ a H , is taken into account. In terms of 
SU(4)xSU(2) L xSU(2) R , H and A acquire VEVs in components of H (1, 1, 1) and ,4(15,1,1), respec- 
tively. Therefore two ^ a (l, 2, 2) need to make singlet by themselves in Z^ a A i ^ a H , but anti-symmetric 
contractions respect to SU(2)p, SU(2)l and SU(2)r indices end up zero. 

* 25 SeeE20]and Appendix H 
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this procedure leads to the following Yukawa couplings: 



Y, 



D 



Qi 
Q2 



( \^ bg-af \ 2 1 



U R\ 



d 5 b g -af Q X 5 
ac—b z q I 



ad b bg-af ^3 
ac—b 2 g 



U R2 



-a-f d 5 



ac—b 2 • 



iS \ 5.5 



ad\ b . 



-id 



+ P H fe«)X 



iS\ \4.5 



ac—b 2 g 

a e -iS 1 a n „iS\\2.5 

-^z^e +p H ge jA 



D c 

1J R3 



9 

bg-af \2 



(5.5) 



Y 1 



E° R1 

E R2 

E c 



R3 \ 



ac—b 2 

Id/A 5 



2 ai ' 





L2 



/We t5 A 4 - 5 



' ac—b 2 1 bg—af b \ \ 4 



+ 



bg-af \2 
9 



(5.6) 



As a result of discrete symmetry, Eq. (15.41) leads to Y u ~ (^-) 2 A 6 , where d q is propor- 



tional to a <5s charge! I Also, Yu turns out to be a real parameters, contrary to Yd 
which possesses nonremovable phases that are converted to the KM phase. In Eq. (15.6[) . 
one can see that the determinant vanishes when 0h is set to zero, i.e., Hd ~ Hd[H 5 }. 
Also note that nonremovable phase in Eq. ( 15.5D vanishes when (3h is set to zero. This is 
the reason why we assigned U(1)a and Z 6 charges so that the down- type Higgs doublet 
is composed as H D ~ H d [H*']+\°- 5 L[C*]. 



5.3 Model analysis 



The resultant Yukawa couplings of Eqs. ( 15.41) . ( 15. 51) . (15. 6 j) have rather restricted forms and 
small numbers of (9(1) coefficients. In this subsection, we examine these Yukawa couplings 
and derive some predictions. Also, we discuss the compatibility with observables. 

First of all, in a semianalytical derivation, where we collect leading contributions 
in each element, Eqs. (15. 4 j ) and (15. 5 p lead to the following Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements. 



V, 



leading 
CKM 



1 

■V,*. 



a(bg-af) 2 



(d q +2d 5 )/3 H 



\V{ 2 V 2S 



2{ac-b 2 ) {(bg-af) 2 +(ac-b 2 )g 2 }p H -a 2 d 2 t 
1 

-^ 23 



rA 



\ 



g ac—b 2 \ 2 
a bg-af 

1 



(5.7) 



* 26 We assume that this nature is maintained even if higher dimensional terms, like Z^ a A 3 ^ a H , is taken 
into account. In terms of SU(4) xSU(2)l xSU(2)r, SU(4) indices of ^ a (4, 2, 2) can not be contracted with 
that of ^0(4, 2, 2), because of the SU(2)f symmetry. Then, it is expected that quark components among 
^0(4, 2, 2) and vE , a (4, 2, 2) pick up at least two Qb charges from the contraction with the adjoint fields 
A(15, 1, 1). Strictly speaking, VEVs of 15-plets that are made by the direct product of adjoint fields can 
give contributions which are not proportional to Qb for quark components of (*a(4, 2, 2) x "1^(4, 2, 2))i 5 . 
However, these contributions generally possess another suppression factors that are determined so that 
irreducible representations do not contain a singlet component. 
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Here, V ub vanishes at the order of 0(X 3 ) and appers at the next order in A in this model. 
It is interesting to note that the current global fit for the CKM matrix elements at low- 
energy [33] or the high-energy extrapolation [25] indeed suggest this type of hierarchical 
structure. 

Also, in the leading order analysis, we found that Eqs.f l5.4p and ( 15. 5 p lead to the 
following relation, 

\V cb Y b \ = \Y C \ (5.8) 

where Y b and Y c denote bottom and charm quark Yukawa couplings. We can guess the 
origin of this relation as follows. The (23) element of Vckm is given as V23 ~ ~ ^33 ' 
but now, Y D2 3 ~ Y U2 2 and Y D33 ~ Y U32 because 5 2 turns into 5 3 , after 5 mixing. Then 
V 23 Y b ~ Y D23 - ^Y D33 ~ y u2 2 y U3 3-y U23 y U32 and thig ig roughly Yc . For comparison with 

the observables, we quote the results of Ref . [25] . 

\V cb Y b \ ~ ^|y c | (5.9) 
6 

where low-energy inputs of V cb , Y b and Y c are extrapolated to the GUT scale. Since 



27 



Eq.f l5.9p is derived in the case of tan /3=10, Eq.f l5.8p fixes tan /3 ~ 6 in the current mode 

Finally, as an example, we numerically calculate the eigenvalues of quark and lepton 
Yukawa couplings, a mixing matrix, and a Jarlskog invariant (Jcp) ; substituting a set 
of trial values for the 0(1) coefficients that appeared in Sec. 5. 2. We set a = 0.6, b = 
_0.5, c = -.7, d 5 = -0.9, d q = 0.4, di = -0.5, / = 1.5, g = -0.9, f3 H = 0.9, 5 = 1.4, 
and repeat the processes described in Sec. 5. 2. We obtain the following eigenvalues of the 
quark and lepton Yukawa couplings, a mixing matrix, and a Jarlskog invariant. 

Y t = 6(5) x 10" 1 Y b = 2(3) x 10" 2 Y T = 3(4) x 10" 2 

Y c = 3(l)xl0" 3 y s = 5(6)xl0" 4 Y^ = 1(3) x 10" 3 (5.10) 

Y u = 4(3) x 10" 6 Y d = 8(3) x 10" 5 Y e = 3(1) x 10" 5 

/ 1 2(2) x 10" 1 2(4) x 10" 3 \ 

\V C km\ = 2(2) x 10- 1 1 10(4) x 10" 2 (5.11) 

\ 20(7) x 10~ 3 10(4) x 10" 2 1 / 



'CP 



1(3) x 10" 5 (5.12) 



Here the parenthetic digits are the corresponding values of Ref. [23] that are extrapolated 
to the GUT scale from the low-energy inputs using the MSSM two-loop RGEfQ From 
Eqs. fl5.10p - fl5.12p . we see that the current model has the capability to reproduce the quark 
and lepton mass spectra, mixing angles, and a Jarlskog invariant. 



* 27 Since uses the MSSM renormalization group equations (RGE), the discussion can be slightly 
changed once heavy fields contributions, that are relevant near the GUT scale, are taken into account. 

*28 p or .j-jjg eigenvalues of down-quark and charged lepton Yukawa couplings, that are originally calcu- 
lated in tan /3 =10, we naively multiply 6/10 to assess the corresponding values in case of tan/3=6. 
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6 Summary and discussion 



In this paper we discussed a model of spontaneous CP violation in the Eg SUSY GUT 
with SU(2)p flavor and anomalous U(1)a symmetries. 

We made a model where CP symmetry is spontaneously broken in the flavor sector in 
order to provide the origin of the KM phase and to evade the SUSY CP problem. However, 
as we saw in Sec J4.2[ a naive construction of the model generally leads to an unwanted 
outcome, Arg[/i&*] = 0(1), when the CP violating effect is taken into account. Then in 
Sec J4.3j we introduced a discrete symmetry in order to cure this difficulty Interestingly, 
this discrete symmetry plays additional roles. It ensures that realistic up-quark mass and 
Cabibbo angle are simultaneously realized without cancellation between 0(1) coefficients. 
Also, severe constraints from the chromo-electric dipole moment of quark can be satisfied 
without destabilizing the weak scale. 

The discrete symmetry reduces the number of free parameters, but the model is capable 
of reproducing quark and lepton mass spectra, mixing angles, and a Jarlskog invariant, 
as we see in Eqs. fl5.10p -f l5.12p . We also obtain the characteristic predictions V u b ~ C(A 4 ) 
and |Vd,Yb| = \Y C \ at the GUT scale. Note that, for the trial set of 0(1) coefficients that 
is used to calculate Eqs.f l5.10~]) - fl5.12p . ArgfDet [Y^/Yd]] is 0(1) and it leads to the strong 
CP phase. However, this phase can be set to zero by the VEV of the axion that originates 
from the U(1)a breaking. 

A more precise fitting of the 0(1) coefficients is beyond the scope of the current 
study, since it requires detailed study of RGEs and threshold corrections. Note that, 
in the model, there are many superheavy fields which affect the RGEs. In particular, 
among the Higgs fields listed in Table [fl positively charged fields C , C, and A' provide 
relatively light fields and these fields start to influence the flow of the gauge couplings 
from the intermediate scales, as depicted in Fig.l of the third paper of Ref. [7]. Then, this 
also results in a change of the flows of the Yukawa couplings. 

In SecJH we discussed the Yukawa couplings of quarks and charged leptons. As for the 
neutrino sector, the left-handed neutrino masses are induced via the see-saw mechanism 
|34j . Here, the right-handed neutrino mass terms are obtained by the following inter- 
actions, ^f^fH^H™, yf^fH^C 1 * 7 , ^f^fC^C™, adding appropriate fields like 
F a , F a , C 27 and C 27 to make these terms singlets under the symmetries listed in Table [TJ 
Then, when we follow the procedure of Ref. [7J , the solar and atmospheric neutrino mass 
squared differences are roug hly derived as Am 2 = C(10" 6 )eV 2 and Am^ tm = C(10" 3 )eV 2 , 
respectively. Note that actual values can change by an order of magnitude, since several 
new 0(1) coefficients appear in the neutrino sector. Therefore, AttIq and Am 2 tm can 
be consistent with the current experimental constraints [33J. This is the same for the 
neutrino mixing angles. In the current model, the Maki-Nakagawa-Sakata (MNS) matrix 
elements [35] are roughly given as 

/ A - 5 1 A \ 
\Vmns\ ~ 1 A ' 5 A - 5 , (6.1) 
\ A - 5 A 1 / 

and these values may vary by an order of magnitude when we replace 0(1) coefficients with 
real numbers. Also, for example, when we change the U(1)a charges as \{c — c— (h — h)} = 
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—0.5, |(Vjvnvs)i3| changes to a more realistic value of C(A 1,5 ). Finally, it is interesting 
to note that, in the model described in Secj5j one of the left-handed neutrino becomes 
massless. 
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A Higgs sector 

In this section, we briefly outline the Higgs sector of the model that is made to leave 
H U [10 H ] and H d [10 H ] + A a5 L[16c] massless. These components are identified as the 
MSSM Higgs doublets Hu and H D , respectively. As stressed in Sec. 5. 2, H D ~ H d [10 H ] + 
A a5 L[16c] is crucial for the nonvanishing KM phase and electron mass. 

In order to distinguish the pair of low-energy massless modes that originate from the 
Higgs sector, here we compose the operator matrix O which induces mass terms among 
the SU(2)l doublet components of H, H, C, C, C, C', A and A] 29 l. O consists of 16 rows 
and 16 columns. In the following, we take a base where indices are in the order 10h, 
10 c , 16c, 16a, 10c, 10c, 16c, 16a', 10s, 10 c , 16 H for rows and 10^, 10 c , 16 c , 16a, 
10c, 10(5' j 16c*" 16a', 10//, 10 c , 16// for columns, respectively. Here we divide O into 
the following 9 blocks: 



10 H 

10f 7 o 



O 



16# V 



4x4 



B 



4x4 



E 



3x4 



.4 



4x4 



a 



4x4 



3x4 



16 



H 







4x3 



D 



4x3 



G 



3x3 



(A.l) 



/ 



where subscripts indicate the number of rows and columns of the partial matrices. Here 
the 0's have no entries and E has a few entries since the corresponding fields are mostly 
negatively charged under U(1)a (SUSY-zero mechanism). On the contrary, the elements 
of C,D,F and G are almost filled since the corresponding fields are mostly positively 
charged. Therefore, in the following, we concentrate on A and B which determine the 
composition of the massless modes. Here B is responsible for the composition of Hu and 
A is responsible for the composition of Hp. 



See 



for more detailed discussion. 



16 



Under the charge assignment given in TableJIJ A and B are filled with the following 
terms. 



A 



10 H 
10c 
16c 
16a 



10c 

( H 2 AC 







10 





c> 



16 



a 



HHC'C 
C'{AFF + Z 4 )C 

C'(AFF + Z 4 )C 

C'AC CAE 



16 A ' 
\ 



A' A ) 



(A.2) 



B 



10 H 10c 16 c 16a 

v(H 2 AC C'H(A + Z 3 )CC CAC'Z 3 \ 

C'{AFF + Zi)C C"(,4 + Z 3 )C 2 C'ACifZs 

HAC 

A' A 



10c 
16c 
16a' \ 



(A.3) 



Then, after symmetry breaking, Eqs.( lA.2l) and (IA.3I) induce the following mass terms 
among SU(2)l components: 
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A 5 


A 13 
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16c 
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A 8 - 5 




16a 
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^ 








A 3 


/ 



(A.4) 



(A.5) 



Note that each (1,1) entry vanishes since A acquires B — L conserving VEV. Therefore, 
Eq. f[Q|) and flAT5|) appropriately realize H D ~ H d [10 H ] + A a5 L[16c] and H v ~ H U [10 H ]. 
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